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Deep inelastic scattering at small x can be described very effectively using saturation 
inspired dipole models. We investigate whether such models are compatible with the 
numerical solutions of the Balitsky-Kovchegov (BK) equation which is expected to 
describe the nonlinear evolution in x of the dipole cross section. We find that the BK 
equation yields results that are qualitatively different from those of phenomenological 
studies. Geometric scaling is recovered only towards asymptotic rapidities. In this 
limit the value of the anomalous dimension y(r,x) at the saturation scale approaches 
approximately 0.44, in contrast to the value 0.63 commonly used in the models. 

At small x, deep inelastic scattering (DIS) can be described as the scattering of a color 
dipole, which the photon fluctuates into, off the proton [2]. The linear BFKL equation, which 
describes the dipole-proton interaction in terms of gluon ladders, predicts an exponential 
growth of the corresponding cross section as log \jx increases, potentially violating unitarity. 
Hence, interactions between BFKL gluon ladders may become important, which leads to a 
nonlinear evolution approximately described by the Balitsky-Kovchegov (BK) equation [3]. 
As a consequence of the nonlinearity, the dipole cross section saturates with decreasing x, 
thereby offering a resolution to the unitarity problem. The inclusive HERA data at low x 
(x < 0.01) could be described well by a dipole cross section of the form a = coA'cbw^ 
where the scattering amplitude A^qbw is given by [4] 

N GBW (r,x) = l-e X p[~\r 2 Q 2 s (x)] , (1) 

r denotes the transverse size of the dipole, <7q ~ 23 mb and the ^-dependence of the satu- 
ration scale is given by Q s (x) — 1 GeV (xq/x) x ^ 2 , where xq ~ 3 x 10 -4 and A ~ 0.3. The 
scattering amplitude depends on x and r through the combination r 2 Q 2 {x) only, which is 
known as geometric scaling and leads to the prediction that the structure function F2 is a 
function of Q 2 /Q 2 {x) only. This prediction was checked in a model independent way [5] and 
holds widely even though the GBW model (fT]) is not applicable at large Q 2 . It should be men- 
tioned that the leading order BK equation leads to a faster evolution in x [6] (Q 2 (x) ~ l/x x 
where A ~ 0.9) than the experimental data seem to favor (A ~ 0.3). This discrepancy can 
be reduced by introducing a running coupling constant. 

Hadron production in d-Au collisions can also be described by saturation inspired dipole 
models [7-9]. However, these data seem to require geometric scaling violation. The dipole 
scattering amplitude modified in this respect is given by [7-9] 

N(r, x) = 1 - cxp \-±(r 2 Q 2 (x)y( r - x A . ( 2 ) 

The exponent 7 is usually referred to as the "anomalous dimension" , although the connection 
of N with the gluon distribution may not be clear for all cases considered below. Following 
partly [7, 10], in [8,9] a few requirements were used to determine a parameterization of 7. 
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Firstly, one assumes that j(r, x) approaches 1 in the limit r — > 0. Therefore the "DGLAP" 
limit N ~ r 2 is recovered for all cc. Secondly, at the saturation scale, r = 1/Q S , 7 should 
be constant to ensure geometric scaling in this region. This constant 7 S is chosen to be 
~ 0.628. The value of 7 S is motivated by a saddle point analysis of a solution of the BFKL 
equation with saturation boundary conditions [6] and also shows up in the traveling wave 
approach [11]. Thirdly, if one writes 7 = 7 S + A7, then A7 should decrease as 1/y for y — > 00 
at fixed r 2 Q 2 . This ensures that geometric scaling is asymptotically recovered. Furthermore, 
the parameters were adjusted in such a way that geometric scaling holds approximately 
for finite y in a growing region between Q s (y) and roughly Qg (y)/ Aqcd • Note that the 
parameterization in [9] is intended to describe N(r,x) in this so-called extended geometric 
scaling region only. To simplify the procedure of the required Fourier transformation of N 
@, 7(r, x) was replaced in [7-9] by "f(l/k,x) where k is the transverse momentum of the 
scattered parton that will fragment into the final state hadron. 

We want to check whether these requirements for 7(7*, x) are compatible with the non- 
linear evolution of the dipole scattering amplitude N. The BK equation for N reads [3] 



dN(r = \x t -y t \,x) _a ± C ^ {xt-yt) 



N{\x t - z t \,x)+N{\z t -y t \,x) 



dy 2it J (x t - z t ) 2 (y t - z t y 

-N(\x t -y t \,x)-N(\St - z t \,x)N(\z t ~y t \,x)y (3) 

Here a s = a s N c /w. We will not consider the impact parameter dependence of N. 

The BKsolver program [12] provides a numerical solution of the amplitude 7V(fc, x) in 
momentum space. In order to use this solution of the BK equation ([3]) to constrain 7(7", x), 
one first has to find N(r,x) by Fourier transforming to coordinate space: 



N(r,x)=r 2 J —± e - tkt ^N{k,x)=r 2 J dkk J (kr)Af(k,x) . (4) 

Using the Ansatz one can extract 7(7*, x) from the resulting N(r,x), 

7 (r, x) = log[log[(l - N(r, x))" 4 ]]/ log[r 2 Q 2 s (x)} . (5) 

This equation requires as a separate input the value of Q s ( x) , which can be found by equating 
the right hand sides of Eqs. ^ and (0| for r = 1/Q S - Combining the resulting values of Q s 
with Eq. (0), we obtain a numerical result for 7(7", x), which is shown in Fig. [1^,. 
The resulting 7(7", x) has the following features: 

1. For r — > 0, 7(7", x) asymptotically approaches 1. 

2. At the saturation scale, 7(7", x) is not a constant. 

3. For decreasing x, 7(7*, x) approaches a limiting curve, r y ao {rQ s {x))^ indicated in Fig.[T] 
by y = 00. Hence, after a longer evolution one indeed recovers geometric scaling. 

The fact that for small distances 7 asymptotically approaches 1 is understandable from 
the BK equation, since in this limit it reduces to the BFKL equation. In the limit of small 
distances, the solution to the BFKL equation is dominated by either the saddle point or 
the initial condition, both leading to 7 — > 1, since here we use the MV model as the initial 
condition, see [13] for details. 
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Figure 1: a) 7(7*, x) resulting from the relations (J4j) and ([5]) as a function of l/(rQ s (x)) and 
y = logxo/x. b) "f(k,x) as a function of k/Q s {x) for various rapidities y = logx^/x. 



It turns out that 7 is clearly not constant, not even at r = 1/Q S , unlike in [7-9]. However, 
asymptotically geometric scaling is recovered as 7 approaches 700. Writing 



j(r,x) = ~foo{rQ s {x)) + Aj(r,x) . 



(6) 



it turns out that, similar to the parameterizations used in [7-9], A-y(r,x) decreases as l/y 
for y — > 00 and fixed r<5 s (x). At the saturation scale 7 is given in the small- x limit by 



lim 7(r 



l/g,(a:),a:)=7 oo (l)w0.44 J 



(7) 



which is significantly below 7 S = 0.628. This is not in disagreement with theoretical expec- 
tations [6, 11]. Rather it indicates that requiring 7 in Eq. @ to be constant at Q s does not 
follow from the BK equation. 

In [7-9], N(r,x) j2]) was considered with 7(7", x) replaced by j(l/k,x). This approxi- 
mation scheme we will discuss next. The procedure of extracting 7 becomes quite different 
when 7 depends on k, since the dipole cross section N then depends on both r and k, so 
that it is not related to M{k, x) by a straightforward inverse Fourier transform ([¥]) anymore: 



M(k,: 



- Mkr) 
r 



1 — exp 



\{r 2 Q 2 s {x))^ 



(8) 



Instead of by using the inverse Fourier transform, we will extract 7 by numerically solving 
Eq. ||HJ|, imposing the following condition. In order to test the Ansatz in [8,9], we will fix 
7(fc,x) in such a way that it equals the constant 7^ ~ 0.628 at the saturation scale. The in- 
dependence of Q s is determined by explicitly solving Eq. ([5]) for k = Q s and j(Q s , x) = 0.628. 
Now we can extract 7 from relation (JSJ) for any given value of x and k. Fig. [1] shows the 
results for 7(A), x) as a function oik/Q s above Q s , for a broad range of rapidities. For small 
rapidities the resulting 7 looks very similar to the one in [9] (cf. Fig. 4 of Ref. [9]). As 
one can see, for larger y the resulting 7 is not compatible with the parameterization in [9] 
anymore; it first decreases before it rises towards 1 asymptotically. 
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For a discussion of additional important issues like the dependence on initial conditions, 
the x-dependence of Q s in our approach and the running coupling case we refer to [13]. 



Discussion & Conclusions 

The numerical solutions of the BK equation do not display exact geometric scaling, although 
they approach a solution showing such scaling at asymptotic y. Assuming the solutions to 
be of the form where scaling violations are encoded in the "anomalous dimension" 7, 
therefore leads to the conclusion that 7(7*, x) is not a function of rQ s (x) exclusively. In 
particular, it is never simply a constant, not even at the saturation scale (r = l/Q s ). At 
asymptotically large rapidities, 7 reaches a limiting function "foo{fQ s (x))- This function 
is universal for a large range of initial conditions [13]. At the saturation scale, 7^ equals 
approximately 0.44, which is considerably smaller than the corresponding values in the 
phcnomenological models [7,9,10]. For small values of rQ s the limiting function seems to 
reach j s [13], in accordance with the traveling wave results of Rcfs. [11]. 

Performing the replacement of 7(7", x) — > 7(l/fc,a;) does allow one to find a solution for 
which 7(fc = Q s ,x) is kept fixed. The behavior of j(l/k,x) is then for small rapidities 
qualitatively similar to the parameterization in [9]. However, the usually considered choice 
"f(k = Q s , x) = 7 S = 0.628 yields some unwanted features, i.e. A7 being negative in a region 
above the saturation scale and the absence of solutions below the saturation scale, although 
the Ansatz was not intended for that region. Keeping j(k = Q s , x) fixed at a smaller value, 
e.g. at 7oo(t , Qs = 1) ~ 0.44, seems more suitable [13], but it remains to be investigated 
whether such a choice allows for a good fit of all relevant DIS, d-Au and p-p data. 

It would be interesting to consider modifications of phcnomenological models for the 
dipole scattering amplitude that are compatible with both the BK equation and the data. 
Given the fact that the BK evolution does not respect geometric scaling around Q s , phe- 
nomcnological parameterizations that reflect this feature would seem a natural choice. For- 
tunately, the LHC and a possible future electron-ion collider will provide data over a larger 
range of momenta and rapidities, so that one can expect to test the evolution properties of 
the models more accurately. 
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